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Ultra high energy gamma rays produce electron-positron pairs in interactions on 
the geomagnetic field. The pair electrons suffer magnetic bremsstrahlung and the 
energy of the primary gamma ray is shared by a bunch of lower energy secondaries. 
These processes reflect the structure of the geomagnetic field and cause experimen- 
tally observable effects. The study of these effects with future giant air shower arrays 
can identify the nature of the highest energy cosmic rays as either 7-rays or nuclei. 
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I. INTRODUCTION 



Ever since the reports of the detection of two cosmic ray showers of energy well above 
10^° eV [jl],|^ the origin and the nature of such events have been subjects of strong interest 
and intense discussion. It is not only very difficult to extend our understanding of particle 
acceleration to such extraordinarily high energies but the propagation of these particles in 
the microwave background and possibly other universal radiation fields restricts the distance 
to their potential sources to several tens of Mpc. 

Conservatively minded astrophysicists are looking for astrophysical sources which may 
contain the environment necessary for stochastic particle acceleration to energies in excess 
of 10^'' eV. Powerful (FRII) radio galaxies @] have been suggested as possible sources. If this 
suggestion were true, the highest energy cosmic rays (HECR) would be most likely protons, 
refiecting the composition of the matter that is available for injection in the termination 
shocks of FRII jets. Others search for powerful astrophysical sources in the cosmologically 
nearby Universe. HECR then could also be heavier nuclei, for which the acceleration is less 
demanding. The propagation of heavy nuclei on short distances (0(10) Mpc) without huge 
energy loss is possible. 

Some cosmologists relate the origin of HECR to topological defects 0. Topological 
defects (TD) scenarios avoid the problems of particle acceleration since they are based 
on 'top-down' evolution. Very massive (10^^ — 10^^ eV) X-particles are emitted by the 
topological defects that later decay into baryons and mesons of lower energy. Most of the 
energy is eventually carried by 7-rays and neutrinos, that are products of meson decay. 
Detected HECR would then most likely be 7-rays. 

Most radically, the origin of HECR has been related to those of gamma ray bursts , 
replacing two extremely luminous mysteries with a single one. In such scenarios HECR are 
most likely to be again protons. We may not be able to observe the sources of HECR since 
every source might only emit a single observed ultrahigh energy particle. 

The nature, the type of the particle that interacted in the atmosphere to generate these 
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giant air showers, could be the key to the understanding the origin of the highest energy 
cosmic rays. The current experimental evidence on the nature of HECR is not conclusive. 
The Fly's Eye experiment, for example, has reported correlated changes in the spectra and 
the composition of the ultra high energy cosmic rays The analysis of the Fly's Eye 



experimental statistics suggests that a change of the chemical composition of the cosmic rays 
from heavy nuclei to protons at ~ 3 x 10^^ eV is accompanied by a change of the spectral 
index of the cosmic ray energy spectrum. One may then conclude that the HECR are 



protons. The other currently running air shower experiment, AGASA, does not observe ||Tl 



such a correlation. A re-analysis of the archival data from the SUGAR experiment [12 
makes the opposite conclusion - a large fraction of the highest energy showers seem to be 
generated by heavy nuclei.. 

A correlation between the arrival directions of HECR with energy > 4 x 10^^ eV with 
the supergalactic plane, that contains most of the galaxies of redshift < 0.03, has been 
reported ||13[. The AGASA experiment has also observed similar correlation in their data 
set , although not fully consistent with the conclusions of Ref . . On the other hand 



the Fly's Eye experiment does not see such a correlation (P. Sommers for the Fly's Eye 
group, private communication). It also has not been observed in the SUGAR data []T3[. 



Even if confirmed in the future, a correlation with the structure of the local universe would 
not answer the question of the nature of HECR. If topological defects are seeds for galaxy 
formation most powerful galaxies and TD would have similar distribution and TD and 
astrophysical scenarios of the origin of HECR are indistinguishable. 

The profile of the 3x10^° eV shower detected by the Fly's Eye develops higher in the 
atmosphere than expected for either proton or 7-ray showers of that energy The 



highest energy shower seen by the AGASA experiment (2 x 10^° eV) exhibits, apart from its 
energy, features that are typical for most of the high energy showers. The currently existing 
air shower arrays cannot drastically increase the experimental statistics and the hope for 
answering the important questions for the nature and origin of HECR is in the construction 
of much bigger shower arrays, such as the Auger project [T^ . 



Even with Auger, however, the nature of HECR will be difficult to study. Shower param- 
eters are subject of strong intrinsic fluctuations and the cross sections that govern inelastic 
interactions at ^/s = 100 TeV are not well enough known. At lower energy (10^^ — 10^^ 
eV) showers generated by heavy nuclei, protons and 7-rays could be at least statistically 
distinguished by their muon content. 7-ray showers have on the average ~ 3% of the muon 



content of proton showers of the same energy |T8| . At ultrahigh energies such approach may 
not be possible - calculations of the muon content of the 7-ray induced showers predict that 
the fraction of GeV muons could be even higher than in proton generated showers [p!9| ,p0(]. 

We suggest a different approach to the study of the nature of the cosmic rays with energy 
above 10^^ eV - to prove (or disprove) that HECR are 7-rays by observing their interac- 
tions with the geomagnetic field. While protons and heavier nuclei are not affected by the 
geomagnetic field, ultra high energy 7-rays interact on it to produce e+e^ pairs. The elec- 
trons themselves quickly lose their energy through magnetic bremsstrahlung (synchrotron 
radiation) before they enter the atmosphere of the earth. Air showers are thus replaced by 
'magnetic + atmospheric' showers that start far away from the surface of the earth and are 
absorbed faster compared to usual air showers. With high experimental statistics one can 
observe the interactions of ultra high energy 7-rays with the geomagnetic field by a study 
of the shower arrival direction in geographical coordinates. If the detected showers do not 
show signs of interactions with the geomagnetic field, the suggestions for 7-ray nature of 
HECR could be proven wrong. 

This article is organized in the following way. Section 2. gives a brief discussion of the 
photon and electron interactions on magnetic fields and of the structure of the geomagnetic 
field. Section 3. describes a calculation of the 'geomagnetic + atmospheric' cascades and 
gives some general results of that calculation. Section 4. calculates shower parameters that 
could be used to confirm the 7-ray origin of HECR and Section 5. contains the conclusions 
from this research. 
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II. PHOTON AND ELECTRON INTERACTIONS IN THE GEOMAGNETIC 

FIELD 

Interactions of photons, and especially of electrons, on magnetic fields have been ex- 
haustively studied because of all the problems they create in particle accelerators. The 
theoretical and some experimental knowledge is reviewed by T. Erber in Ref. 

Magnetic pair production is guided by the parameter = [l/2][hv/ mc^] [B± / Bcr] , where 
Bcr = m^c^ /eh = 4.414 x 10^^ Gauss and B± is the component of the magnetic field that 
is normal to the 7-ray trajectory. The 7-ray attenuation coefficient, i.e. the fraction of 
photons that undergo pair production in magnetic field of strength B± per unit distance is 
given by 

a,(T,) = 0.16^^ir,%(2T,/3)cm-^ (1) 

The maximum attenuation is reached at 7-ray energy of 121710^ {B / Bj_) while the cross 
section of the process is linearly proportional to the magnetic field strength B±. 

Similarly the magnetic bremsstrahlung (synchrotron radiation) is guided by Tg = 
[E/mc^][B±/ Bcr]- The radiation emitted by an electron of energy Ee in magnetic field 
B± per unit distance is distributed as 

where J = [hu/Ee][l + hu/Ee]/3Te. 

To demonstrate the strength of the 7-ray interactions in the geomagnetic field we show 
in Fig. ID the distributions of the distances from the surface of the earth at which 7-rays of 
different energy pair produce. The 7-ray trajectory is taken to be normal to the field lines 
of a magnetic dipole centered at the center of the Earth with magnetic moment of 8.1 x 10^^ 
Gauss/m. One could see that the 7-rays of the energies of interest interact in a relatively 
narrow range of distances not further than 3R^. The narrow peak plotted at altitude of 20 
km represents 7-rays that survive, i.e. interact in the atmosphere before they interact in 



5 



the geomagnetic field. 12% of tlie 7-rays witli energy 10^° eV (and none at liiglier energy) 
survive. 

Tfie spectra of tlie 7-rays emitted in magnetic bremsstralilung depends quite strongly 
on the magnetic field strength. For strong fields the energy distribution of the secondary 
photons is quite fiat. Fig. ^ shows the energy loss of 10^° electrons in magnetic fields of 
strength log^^Q B± = -0.5, -1, -1.5, etc. Gauss as a function of the secondary photon energy. 
In the dipole field model described above a field of 0.1 Gauss corresponds to a distance of 
0.468 i?0 above the surface of the earth, and 0.032 Gauss - to 1.15 R^. These distances cover 
much of the primary 7-rays interaction range shown in Fig. [l|. Since lower energy 7-rays 
pair produce close to the earth, the magnetic bremsstrahlung of their secondary electrons 
is harder. The energy spectra of the 7-rays in the bunch that enters the atmosphere after 
pair production and magnetic bremsstrahlung tend to be almost independent of the primary 
7-ray energy E^. 

A. Structure of the geomagnetic field 

Gamma rays arriving at any experimental location under different zenith (^9) and az- 
imuthal (0) angle will see a different geomagnetic field. They will thus cascade differently 
before reaching the atmosphere. At small i), close to the vertical direction, the variation 
with (f) is insignificant. At relatively large i), more than 30°, the field strength for most 
locations changes by factors of 3 or more for different values of (p. 

A more quantitative calculation of the strength of the field encountered by the incoming 
7-ray is trivial for any model of the geomagnetic field but has to be performed for each 
location, and separately. We have attempted to obtain a slightly more general result 
for several experimental locations. Fig. |^ shows the transverse component {B±) of the 
geomagnetic field as a function of the azimuthal angle at which it arrives to the detector. 
Since the variations for different zenith angles have the same aspect, we have integrated 
over ?9 from to 60°, weighting the field values with the solid angle. The 1991 International 
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Geomagnetic Reference Field model (IGRF) is used for this calculation. 

Four of the locations for which B± is shown are in the Northern hemisphere and only 
one (Sydney, shown with dash-dash) is South of the equator. Since the smallest B± is seen 
in the direction of the magnetic pole that is closer to each location, northern and southern 
locations have opposite field strength dependence on cj). At the moderate latitudes of these 
detector locations, the detailed differences between the Northern hemisphere detectors are 
minor. (A detector located at the geomagnetic equator would have a symmetric response 
to geomagnetic North and South directions). The difference between the maximum and 
minimum field strengths is almost a factor of 5. 

For each one of these detectors, as well as for any other detector location, one could 
determine a region in azimuth, where the field strength is the lowest and the incoming 7- 
rays would be affected minimally by the geomagnetic field and a region where the effect 
of the geomagnetic field is at maximum. For the location of Sydney, e.g. 7-rays arriving 
with 130° < (f) < 215° would see B^^ <0.02 Gauss at a distance of 1 i?® and 7-rays with 
255° < 4> < 90° would see more than 0.04 Gauss at the same distance. The idea is that 
7-ray fluxes arriving from these two regions may have observable characteristics that are 
different enough to be distinguished experimentally. We continue to study the cascading of 
ultra high energy gamma rays in geomagnetic fields with different strength, corresponding 
to these two regions. 

III. CASCADING IN THE GEOMAGNETIC FIELD 

We simulate the electromagnetic cascading in the geomagnetic field by injecting 7-rays 
of energy at a distance of 5 i?© from the surface of the earth on a trajectory with 
angle i? relative to the vertical direction at the intersection with the surface. The 7-ray is 
propagated with a stepsize Ax (from 1 to 10 km) until the 7-ray pair produces or reaches 
the atmosphere. The atmosphere is defined to be at altitude of 20 km above the earth's 
surface. Gamma rays that reach the atmosphere 'survive' and interact in the atmosphere to 
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produce air showers with their original injection energy. 

If the 7-ray produces an electron-positron pair, the pair electrons are followed in similar 
way, by calculating their radiation spectrum on every step of propagation. The synchrotron 
7-rays are tabulated in energy, starting at 10^^ eV. The assumption here is that secondary 
7-rays of energy less than 10^^ eV do not contribute significantly to the cascades that are 
observed deep in the atmosphere. This lower energy end of the magnetic bremsstrahlung 
spectrum, as well as the electrons of energy below 10^"^ eV that enter the atmosphere, contain 
always less than 2% of the primary 7-ray. 

Each particle produced in the geomagnetic field, as well as the 'surviving' primary 7-rays 
then generate atmospheric cascades. The profiles of these cascades are added up to calculate 
the composite shower profile, generated in the atmosphere by the injected primary 7-ray or 
the products of its interaction in the geomagnetic field. 

The actual calculation is performed using the dipole magnetic field model with magnetic 
moment of 8.1 x 10^^ Gauss/m with two scale factors of 0.25 (low field) and 1.25 (high field), 
which generate magnetic field strengths approximately equal to the maximum and minimum 
values shown in Fig. ^. 

To study the 'survival' probability in these two field models we made calculations for 
two extreme zenith angles: = 0° and 60°, which is the maximum zenith angle at which 
air showers can be reliably detected and analyzed. Fig. ^ shows the survival probabilities 
at high (lefthand strip) and low magnetic field strengths. The lefthand boundary of each 
strip corresponds to propagation at ^9 = 60° and the righthand boundary is for = 0°. 
7-rays approaching the earth at higher zenith angles spend significantly more time in higher 
geomagnetic field strengths and have a higher interaction probability. The lefthand edge 
of the high field strip and the righthand edge of the low field strip practically bracket the 
survival probability space for 7-rays approaching any location at the earth surface with 
zenith angles smaller than 60°. 7-rays arriving at higher angles may be absorbed faster. 

Several calculations of the 7-ray cascading in the geomagnetic field have been previously 
performed [^,^,^]. Our results are in a good agreement with the main results of all 
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of them. Our calculation is generally a refinement of previous ones, which nevertheless 
reveals some practically important features in the cascading process. Previous calculations 
conclude that there will be a 'cutoff' in the energy spectrum of the 7-rays that reach the 
atmosphere, because of the very soft spectrum of the secondary photons, generated by 
magnetic bremsstrahlung. This conclusion is partially due to the relatively rough treatment 
and low statistics in the previous work. Fig. |^ shows the number of secondary 7-rays and the 
energy that they carry. E'^ = 3 x 10^° eV in this example. Although the number of secondary 
7-rays of energy above 10^^ eV is on the average only 6.3, they carry 48% of the primary 
energy. This is also important for the development of the subsequent air showers, because 
at energies above 10^^ eV the LPM effect ||2^, which suppresses the electromagnetic cross 



sections at high energy and slows the development of the air showers, becomes important in 



air. 



IV. ATMOSPHERIC SHOWERS 

Gamma rays of energy above 10^^ eV, if they do exist, would only be detectable by giant 
air shower arrays located on the surface of the earth. Air shower arrays consist of a large 
number of counters that trigger in coincidence when the shower from arrives. The shower 
direction is determined by the arrival time of the shower front at the different counters. A 
fit of the density in the separate counters reconstructs the total number of shower particles, 
the showers size N^, which is then used to determine the primary energy. 

The output of our Monte Carlo simulation includes the shower sizes calculated for several 
atmospheric depths from the cascading of all secondary (and primary, if the injected gamma 
rays did not pair produce) 7-rays in the atmosphere. The profiles from individual secondary 
7-rays of energy above 10^^ eV are calculated with an account for the LPM effect, although 
the effect is not significant below 10^^ eV. The depths are arbitrarily chosen to include a 
realistic range for typical large air shower experiment and correspond to an array at vertical 
depth of 860 g/cm^ and zenith angles with cos{} = 0.9, 0.8, 0.7, 0.6 and 0.5. Because the 
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development of purely electromagnetic showers in the atmosphere (apart from their muon 
content) does not depend significantly on the atmospheric density profile, the examples given 
below could be scaled also to different altitudes and zenith angles. 

Fig. 1^ shows a general and important shower parameter - the average size (< A'^e >) 
generated by 7-rays of different energy. The solid lines are for low field (scale factor of 0.25) 
and the dashed lines are for high field (scale factor of 1.25). From top to bottom the lines 
show Ne{E^) at the 5 different depths of 956, 1075, 1229, 1433 and 1720 g/cm^. Except for 
the deepest observation level, Ng is multiplied by the factor shown by each curve to make the 
figure readable. In the absence of interactions in the geomagnetic field, and for lower E^, the 
shower size have a power low dependence on E^, = E^ with a >1. The power law index 
a depends on the column density between the depth of the shower maximum X^ax and the 
detector. The size at maximum N^ax is exactly proportional to E^ and a is bigger than 
unity because the depth of shower maximum grows with energy as X^ax = log(£'°/81MeV) 
radiation lengths (1 r.l. = 37.1 g/cm^ in air). 

The dependence shown in Fig. ^ is more complicated because in this energy range showers 
are already at or before their maximum development at some of the shallower observation 
levels. The role of the magnetic field strength on the Ng dependence on is easier to 
understand for the deepest levels of observation. Compare, for example, the two curves for 
depth of 1720 g/cm^ with the 7-ray survival probability of Fig. ^ At low energy, where 
there are no interactions on the geomagnetic field, the two curves are the same. The solid 
curve (low field) starts bending at 7-ray energy 2 x 10^*^ eV where the primary 7-rays start 
interacting in the geomagnetic field. Because of these interactions the primary 7-ray is 
replaced by a bunch of 7-rays of lower energy. The composite shower reaches maximum 
at shallower atmospheric depths and is significantly absorbed at the deep observation level. 
The same happens at energy lower by about one order of magnitude in the high field case. 
Although it is outside of the energy range of Fig. ^, at some higher energy, where all 7-rays 
interact on the geomagnetic field, the two curves will join again. 

To explain the behaviour at the shallow observation levels one has to take into account 
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some of the details of the cascading in the geomagnetic field, namely the shape of the energy 
spectra of the secondary photons as a function of field strength, which is shown in Fig. |^. 
Although the primary 7-rays interact in the same way, in the high field case the energy 
spectra of the secondary 7-rays are harder, hard enough to generate showers that are not 
absorbed at the level of 956 g/cm^. One could hardly see a tiny deviation of the strong 
field (dashed) curve in the region of -E^ = 3 x 10^^ eV. At higher energies the secondaries 
are energetic enough to produce dependence very close to a power law. When the 
primary 7-rays start interacting in the low field, however, the picture is slightly different. 
The secondary 7-ray spectra are softer, the composite showers reach maxima at shallower 
depths and are correspondingly absorbed when they reach the observation level. The two 
curves will join asymptotically. 

All other levels show intermediate behaviour where the relation between the depth of 
observation and Xmax also contributes to the exact shape of the curve. 

Fig. ^ shows the strong differences in the observable parameter Nf. which is introduced 
by the strength of the geomagnetic field. It can not be used, however, for analysis of 
experimental data because is not a directly measurable parameter. What experiments 
can do, and usually do, is to produce a spectrum of the measured shower sizes N^. Such 
spectra for the three deeper observation levels are shown in Fig. |^. The solid histogram 
corresponds to the low field and the dashed one - to the high field case. 

The histograms are result of a simulation, where E'^ is sampled from a (-E^)"^ differential 
primary spectrum between 10^^ and 10^^ eV. At low N^, the spectra are always higher for 
the low field case, including the two observation levels that are not shown in Fig. ^. At the 
high Ne side and for shallow observation levels the high field case shows higher spectrum, as 
could be expected by the results shown in Fig. ^ and as seen for the shallowest level plotted 
in Fig. 1^. The biggest difference is at the deepest observation level, where the spectra are 
different as much as a factor of 10. The differences between the size spectra decreases for 
shallower observation levels, and is probably not detectable for the two shallowest levels, 
which are not shown. 
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V. DISCUSSION AND CONCLUSIONS 



The size spectra of Fig. |^ show that it is possible to detect the difference between a flux of 
7-rays that reach the earth after cascading in geomagnetic field of different effective strength. 
In practical terms this means that any experiment that is able to collect large enough 
experimental statistics should see different N^. spectra in different azimuthal directions if 
HECR are indeed 7-rays. We have not attempted to look for this effect in the existing 
experimental statistics, because it is not large enough to reveal such effects. 

The Auger project is an entirely different story. It proposes the construction of two 
air shower arrays, at least 3 000 km^ each, in the Northern and in the Southern hemisphere. 
For comparison, the area of the largest current detector (AG AS A) is 100 km^. Each one of 
these detectors will have the collecting power of more than 5 000 showers above 10^^ eV per 
year. If the HECR primaries are 7-rays the Norther hemisphere detector should see spectra 
similar to the low field case of Fig. |^ in northern direction (the exact direction and value of 
the minimum field depends on the location of the array) and the strong field case in showers 
coming from South. The southern detector will have the opposite effect. If this were the 
case, the difference of the size spectra would prove that the HECR are gamma rays. 

The actual effects could even be stronger than shown in Fig. ^ because the simulation 
on which it is based propagated all primary 7-rays along the vertical {d = Q°) direction. In 
fact, as shown in Fig. ^, the increase of the interaction probability in the geomagnetic field 
increases by a non-negligible factor when the exact particle trajectory is accounted for. Our 
air shower simulation also does not account for the magnetic bremsstrahlung of the shower 
electrons which at high E^. and low atmospheric density < 10^^ g/cnr" could be important 
and could accelerate the shower development. 

In principle the interplanetary magnetic field has to be added to the 'target' magnetic 
field. Gamma ray arriving from a cone centered on the Sun would be absorbed far away from 
the earth and possibly not detectable. The sun could thus be visibly in ultra high energy 
7-rays. The exact dimensions of the region where 7-rays are absorbed in pair production on 
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the solar magnetic field, carries valuable information on the magnetic field in the vicinity of 
the sun. This is an interesting although purely academic problem, because the statistics of 
such events is always going to be negligible. 

Although we have not done it for this paper, there will be effects, similar to the A'^e 
ones, on the muon content of the 7-ray initiated air showers. It is well known that at these 
extremely high energies the number of soft muons (0.3 to 2 GeV) in 7 initiated showers 
is comparable to this of hadronic showers [^,^. The number of soft muons has an 
dependence very similar to A^g, because the low energy muons decay readily when Xmax is 
distant from the observation level. The decay length of 1 GeV muons is ~6 km. A picture 
similar to the A^g spectra in Fig. |^ will develop as a result of the cascading in geomagnetic 
fields of different effective strength. The major difference between the behaviour of the 
electron size and the muon size is that attenuates as a function of the column density, 
while A^ attenuates as a function of the distance, i.e. A^ will depend strongly on the shower 
zenith angle 

The study of the two major components of the giant air showers can reveal the nature 
of the highest energy cosmic rays. If the specific dependence on the shower arrival direction 
is observed, then the highest energy cosmic rays are 7-rays. A non observation of this effect 
would leave us with the choice between protons and heavy nuclei. 
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FIGURES 

FIG. 1. Distribution of the interaction points of 7-rays of energy 10^^ eV (solid line), 3. 16 x 10^° 
eV (dotted line) and 10^'' eV (dashed line). The interaction points are the vertical distances from 
the surface of the earth. The shading in the left hand side of the figure represents the atmosphere. 



FIG. 2. Energy loss of 10 eV electrons as a function of the strength of the magnetic field 
and the energy of the secondary photons. The field strength is indicated by the respective curve 
as log]^o(-^-L/Gauss). 

FIG. 3. The strength of the geomagnetic field component that is perpendicular to the 7-ray 
trajectory as a function of the azimuthal angle ^ at which the particle arrives at the location is 
shown for a distance of 1 i?© from the detector. The field strength is integrated over zenith angles 
1? from to 60° accounting for the solid angle. The calculation is performed for the locations of 
several air shower arrays: a.) Fly's Eye (40N, 112W) - solid line; b.) Yakutsk (62N, 129E) - dots; 
c.) Akeno (35N, 138E) - dashes; d.) Haverah Park (54N,2W) - dash-dot; e.) Sydney (30S, 150E) 
- dash-dash. The 1991 IGRF model of the geomagnetic field is used in this calculation. 

FIG. 4. Survival probability for 7-rays of energy between 10^^ and 10^^ eV in the dipole 
geomagnetic field model described in the text with scaling factors of 0.25 (righthand strip) and 
1.25 (lefthand strip). The lefthand edge of each strip shows the survival probability for 7-rays 
approaching the surface of the earth with a zenith angle = 60° and the righthand edges are for 
^ = 0°. 

FIG. 5. Energy spectrum of the secondary 7-rays that reach the atmosphere after the cas- 
cading of a primary 7-ray of energy 3 x 10^° eV - solid line, lefthand scale. The dotted histogram 
and the righthand scale show the amount of energy carried by the secondary 7-rays in each bin. 
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FIG. 6. Relation between the average shower size iVg and the primary 7-ray energy for 
the five observation levels defined in the text for cascading in high (solid) and low (dashed line) 
strengths of the geomagnetic field. values are multiplied by the factor indicated by the curves. 

FIG. 7. Integral shower size N^, spectra generated by primary 7-rays sampled on a {E^)~'^ 
differential spectrum between 10^^ and 10^^ eV. The solid histograms show the low field case and 
the dashed histograms are for high field. The observation levels are 1720, 1433, and 1229 g/cm^ 
from left to right. 
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